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We show how, within a preformed pair scenario for the cuprate pseudogap, the nodal and antinodal responses 
in angle resolved photoemission spectroscopy necessarily have very different temperature T dependences. We 
examine the behavior and the contrasting T dependences for a range of temperatures both below and above T c . 
Our calculations are based on a fully microscopic T-matrix approach for addressing pairing correlations in a 
regime where the attraction is stronger than BCS and the coherence length is anomalously short. Previously, the 
distinct nodal and anti-nodal responses have provided strong support for the "two-gap scenario" of the cuprates 
in which the pseudogap competes with superconductivity. Instead, our theory supports a picture in which 
the pseudogap derives from pairing correlations, identifying the two gap components with non-condensed and 
condensed pairs. It leads to reasonably good agreement with a range of different experiments in the moderately 
underdoped regime and we emphasize that here there is no explicit curve fitting. Ours is a microscopic rather 
than a phenomenological theory. We briefly address the more heavily underdoped regime in which the behavior 
is more complex. 

PACS numbers: 03.75.Hh, 03.75.Ss, 74.20.-z 



I. INTRODUCTION 
A. Background Literature 

An important dichotomy is emerging in descriptions of 
the mysterious pseudogap phase of the cuprates which has 
resulted in different theoretical scenarios^. At the heart of 
this dispute is whether the pseudogap observed in the normal 
state is derived from the superconductivity itself or whether 
it results from a competing, but somewhat elusive, order pa- 
rameter. Experiments (i) which directly study this anoma- 
lous normal phase have provided evidence for both points 
of vie w . 2 ' 3 ' 4 ' 5 . However, there is an even larger class of 
recent experiments (ii) which address the superconducting 
phase. These are based on angle resolved photoemission 6,7 ' 8 
and Raman scattering 9 ^ as well as scanning tunneling 
microscop y 11,12,13 ' 14 . They quite generally reveal that there 
are two distinct temperature dependences associated with the 
behavior of the spectral function and related properties, in the 
nodal and antinodal regions of momentum space. The nodal 
response appears to reflect superconducting order whereas the 
anti-nodal response is much less sensitive to T c . For this 
reason, it is speculated, that the pseudogap may derive from 
a competing order parameter. Finally, there is a third class 
of experiments (iii) which probe the behavior as the system 
evolves from above to just below T c and establish that the 
transition is clearly second order. Here, for example, one sees 
a very smooth evolution of the ARPES response in the anti- 
nodal direction 1516 . Many other propertie a 17 ' 18 which depend 
on the excitation gap show no clear signature of T c . This 
is generally interpreted as evidence in favor of a precursor- 
superconductivity origin to the pseudogap. 

It is the last two classes of experiments which are the fo- 



cus of this paper. Indeed, there is very little in the theoretical 
literature which addresses these phenomena. Rather the em- 
phasis has been on the ground state or on the normal, pseu- 
dogap phase. Our goal is to show how to reconcile, in partic- 
ular, the experiments of class (ii) with a preformed pair sce- 
nario. Moreover, it is possible that the arguments presented 
here can be viewed as "modular" in the sense of applying 
to alternate precursor superconductivity approaches such as 
the "phase fluctuation" approach 19 or the RVB scheme 2 ^. We 
stress that there appear to be no counterpart studies of the 
intermediate temperature broken symmetry state within the 
more widely espoused phase fluctuation scheme 19 . Our expla- 
nation of the dichotomy is built around a picture in which the 
short coherence length cuprates are somewhere between BCS 
and Bose-Einstein condensed (BEC) systems. This crossover 
scheme seems to be gaining in support^ 2 !, and is now widely 
studied in the cold Fermi gase s 17 ' 22 ' 23 . Our emphasis here is 
on moderately underdoped cuprates where at the lowest tem- 
peratures the spectral properties appear to conform to that of 
a simple rf-wave BCS-like state 7 -" 2 ^. While the behavior ap- 
pears to be much more complex in the heavily underdoped 
regime, nevertheless, there is a smooth evolution with dop- 
ing and all the indications for distinct nodal and anti-nodal 
responses are present at moderate underdoping. Thus, we feel 
the same qualitative physics regarding the origin of the pseu- 
dogap is appropriate to both moderately and heavily under- 
doped cuprates. 

We build on a d-wave BCS-like ground state where the 
variational parameters are determined in conjunction with a 
self-consistency condition for the chemical potential, fi. This 
self consistent treatment of fi (which is close to but differ- 
ent from Ep) is necessary 25,26 to accommodate the relatively 
short coherence length of the cuprates. Our contribution in 
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the past 17,27 has been to address the associated finite temper- 
ature behavior within a microscopic, diagram-based T-matrix 
theory. In earlier papers the anomalous behavior of the Nernst 
coefficient and of the optical conductivity were also addressed 
within this framewor k 28 ' 29 , along with other experiments 17 , 
including 30 the nature of the specific heat jump and the behav- 
ior of the conductance dl/dV. Moreover, a number of years 
ago^S we presented a description of the spectral function with 
special emphasis on how superconducting coherence would 
be evident in the presence of a normal state pseudogap. A 
central point of the present paper is to show that these calcu- 
lations (which predate the actual experiments 7,24 by five years 
or more), yield very good semi-quantitative agreement with 
a wide range of more recent ARPES experiments without in- 
voking any fitting parameters or phenomenology. 

At the onset, we present the simple physical picture of the 
different ARPES spectral gap responses as a function of k. 
We note that the nodal regions are associated with extended 
gapless states or Fermi arcs 24 which are now rather reasonably 
well understood 31 within a pre-formed pair scenario above T c . 
Their collapse below T c has also been addressed within the 
present formalism 32 . One can anticipate (as we find) that the 
arcs are sensitive to the onset of the order parameter, which 
we call A sc , in the same way that a strict BCS superconduc- 
tor, (which necessarily has a gapless normal state), is acutely 
sensitive to the onset of ordering. By contrast, the anti-nodal 
points are not as affected by passing through T c because they 
already possess a substantial pairing gap in the normal phase. 
One will also reach this conclusion by arguing that it is a 
corollary of a second order transition. If there is a difference 
between the nodal and anti-nodal responses above T c (as is 
implicit in the presence of the Fermi arcs), it must persist, as 
we find here, for some range of temperatures below T c . A key 
point to implementing this physical picture is the realization 
that the excitation gap which we call A is, at all temperatures 
(except strictly T = 0), different from the order parameter 
A sc . This distinction trivially holds in the normal, pseudogap 
phase. 



B. Physical Picture of BCS-BEC Crossover Scenario 

Before presenting our microscopic scheme it is useful to 
sketch a simple physically intuitive approach of the BCS- 
BEC crossover scenario at finite temperatures. This approach 
should be seen to be distinct from the phase fluctuation sce- 
nario. As shown in Figure Q] the precursor superconductiv- 
ity here refers literally to pre-formed pairs, rather than (as 
in the phase fluctuation scheme 19 ) to extended regions or 
grains where the order parameter amplitude is well established 
while the phase is uncorrected. These pre-formed pairs arise 
from a stronger-than-BCS attraction. This strong attraction 
breaks the usual degeneracy between A and A sc or the sim- 
ilar degeneracy between the pair formation temperature T* 
and condensation temperature T c . Within this BCS-BEC sce- 
nario, the mechanism for pairing need not be specified. The 
physics focuses on the anomalously short coherence length 
of the cuprates (associated with strong attraction or high T*), 
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Figure 1: (Color online) Cartoon of the model showing non- 
condensed pairs in red, open ellipses and condensed pairs in blue, 
closed circles. The number of non-condensed pairs scales with the 
height of the red region in the following figure. 
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Figure 2: (Color online) Contrasting behavior of the excitation gap 
A(T) and superfluid order parameter A SC (T) versus temperature. 
The number of noncondensed pairs varies as A^ g = A 2 — Ag C . 



whereas in the phase fluctuation scenario the focus is on the 
anomalously low plasma frequency- leading to soft phase 
fluctuations, and more mesoscopic regions of superconduc- 
tivity. 

Figure Q] shows the schematic behavior as one passes from 
above T* to the fully condensed ground state. The red (dot- 
ted) lines enclose Cooper pairs with net finite momentum, 
while the blue (solid) pairs correspond to the components of 
the condensate which are at zero center of mass momentum 
and have phase coherence. The third panel with < T < T c 
is the most interesting from the perspective of the present pa- 
per. This is the regime about which there has been very little 
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theoretical discussion in the literature and this is the regime 
where the interesting "two-gap" scenario physics is emerg- 
ing. Here one sees a three-way co-existence: of the conden- 
sate, the fermionic excitations (denoted by a single spin arrow) 
and of pair excitations or non-condensed pairs. When there 
is a stronger than BCS attractive interaction, preformed pairs 
above T c , which are responsible for the pseudogap, do not 
disappear, but rather evolve smoothly below T c into this new 
form of condensate excitations arising from non-condensed 
pairs. This leads to two gap contributions 33 in the superfiuid 
phase representing the finite momentum pair excitations of the 
condensate (associated with the component, A pg ) and the con- 
densed pairs (associated with the order parameter, A sc ). 



In this two-gap preformed pair scenario there is a gradual 
inter-conversion of non-condensed to condensed pairs as the 
temperature decreases. This is shown in Figure [2] where the 
energy gap parameters are schematically plotted. Above T c 
but below T* the excitation gap reflects the fact that one has 
to add energy in order to create fermionic excitations or break 
pairs. This excitation gap A smoothly evolves below T c as 
in a second order phase transition, while precisely at T c the 
order parameter A sc opens up. The difference between the 
(squares) of these two parameters can be associated with the 
number of non-condensed pairs. Figure|2]thus shows that the 
number of non-condensed pairs is finite below T c provided 
the temperature is different from zero. We will show, using 
our microscopic scheme that the two gap components add 
in quadrature 33 to yield the thermodynamical gap parameter 
A(T). Importantly, A(T) is essentially temperature indepen- 
dent as a consequence of this inter-conversion from A pg (T) 
to A SC (T). Just as there are two gap parameters, there are 
two temperature scales: T* marking the gradual onset of the 
pseudogap, as well as T c which marks the appearance of the 
condensate. 



How do we understand the phase diagram of the cuprates 
within the BCS-BEC crossover approach? Our interest here 
is not on the details of the hole concentration dependence 
although this has been discussed elsewhere 22,34 . There is 
a pronounced competition between T* and T c within the 
BCS-BEC crossover scenario, as the attractive interaction \U\ 
increases 35 . Indeed, when T* increases (as for example with 
underdoping), T c will ultimately decrease. This is due to the 
fact that at large \U\, it is energetically very expensive to un- 
bind a pair of fermions, as is required in the pair hopping pro- 
cess, a large effective pair mass is then responsible for a small 
T c . In the d-wave case 17 this pair hopping is even more re- 
stricted because of the extended size of the pair, which leads to 
pair localization, and quite possibly the "singlet glass" phase 
which has been reported recently 13 . Importantly, this con- 
comitant cessation of T c occurs while the system is still deep 
in the fermionic regime where the chemical potential fi is pos- 
itive, suggesting a phase diagram not so different from that of 
the cuprate s 22,36 . 



H. OVERVIEW OF FULLY MICROSCOPIC THEORY 

Having discussed the simple physical picture we next re- 
view in more detail the underlying microscopic (T-matrix) 
theoretical formalism, which leads to i t 17 ' 30 i 36 . 

A. T-matrix Theory 

We begin with a BCS-like ground state: v&q = IIk(wk + 
Uk<4 t c _fc 4)|0)' where the parameters Uk and Vu are deter- 
mined variationally in conjunction with a self-consistent con- 
dition for the chemical potential, fi. Knowing, as we now do, 
that at the lowest temperatures the spectral properties appear 
to conform to that of simple BCS-like d-wave pairing serves 
to justify this starting point. We have extensively addressed 
the finite temperature behavior associated with this fully con- 
densed ground state as well as the spectral properties 30 . 

To address d-wave pairing in the cuprates we need to in- 
corporate specific k dependent factors so that the gap param- 
eters in the self energy acquire the form Ak, sc = A sc </?k 
and Ak lPg = A pg Lpk, where we introduce </?k = cos(20), 
to reflect the d-wave k dependence along the Fermi sur- 
face. We adopt a tight binding model for the band disper- 
sion ek = 2£(2 — cosk x — cosky) + 2t z (l ~ cosk z ) + 
4i'(l — cos k x cos k y ). It should be stressed that all gap param- 
eters have the same k dependence. The additional effects of 
anisotropy (beyond those in ip^) which appear in the measured 
spectral gaps, are not presumed to be present in the initial gap 
parameters. 

We will next briefly summarize the key equations which 
emerge from our T matrix scheme 17,27 . Throughout this pa- 
per, we adopt a four- vector notation: Q = (ifi;,q), K = 
(iw n ,k), and J2 Q ee TJ2i E q . Ek = T E„Ek> where 
uj n and SI/ are the odd and even Matsubara frequencies, re- 
spectively. We also take h — ks — 1. Within the present 
approach there are two contributions to the full T-matrix 

t tpg ~\- t SC (1) 

where 

t sc {Q) = -5jr8(Q). (2) 
Similarly, we have two terms for the fermion self energy 

= ]Tt(Q)Go(Q-XVk- q /2, (3) 

Q 

where Go is the bare Green's function. It follows then that 

A 2 sc A 2 sc 

E ac (k,iu; n ) = — ^ = - (4) 

iu n + e k - fi iuj n + 4k 

Here 4k = £k — Throughout, the label pg corresponds to 
the "pseudogap" and the corresponding non-condensed pair 
propagator is given by 
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where the pair susceptibility x(Q) has to be properly chosen 
to arrive at the BCS-Leggett ground state and U is the attrac- 
tive pairing interaction. We impose the natural BEC condition 
that below T c there is a vanishing chemical potential for the 
non-condensed pairs 



Impair 0) 

which means that t pg (Q) diverges at Q = 
Thus, we approximate 33,37 Y, pg (K) to yield 



(6) 

when T < T c . 



E pg (K) a -G (-K)Al 



py 



(T<T C 



(7) 



with 



normal state (where n pa i r is nonzero), Eq. dT) is no longer a 
good approximation, although a natural extension can be read- 
ily written dowiA 

To evaluate A pg in Eq. (0 we note that at small four-vector 
Q, we may expand the inverse of t pg after analytical contin- 
uation. Because we are interested in the moderate and strong 
coupling cases, where the contribution of the quadratic term 
in £1 term is small, we drop this term and thus find the follow- 
ing expression, which, after analytical continuation, yields the 
expansion 



tpgiQ) — 



1 



(15) 



E WQ)- 



(8) where il^ = q 2 /(2M*) and where Z is the inverse residue 
given by 



It follows that we have the usual BCS-like form for the self 
energy 



E(k,iu„) 



with Ak = A(^k and 



A k 



(T < T c ) 



A 2 (T) = A 2 (T) + A 2 C (T). 



(9) 



(10) 



As is consistent with the standard ground state constraints, 
A pg vanishes at T = 0, where all pairs are condensed. 

Using this self energy, one determines G and thereby can 
evaluate t pg . Then the condition that the non-condensed pairs 
have a gapless excitation spectrum (fJ. pa ir = 0) becomes the 
usual BCS gap equation, except that it is the excitation gap A 
and not the order parameter A sc which appears here. We then 
have from Eq. d6) 



1 + 17 



1 



k 



2/(£k) 2 

<Pk = o, 



2£ k 



T<T C , 



(ID 



ylk + A 2 , is the quasiparticle dispersion. 



where 

To close the loop, for consistency we take for the pair sus- 
ceptibility 



x(Q) 



E 

K 



HereG = [G^-Y,] 
using 



G Q (Q - K)G(K) v l_^. (12) 
_1 is the full Green's function. Similarly, 



= 2]TG(iq 



K 



one derives 



E 



+ 2^f(E k ) 

^k ^k 



(13) 



(14) 



which is the natural generalization of the BCS number equa- 
tion. The final set of equations which must be solved is rather 
simple and given by Eqs. ([8]), ( fTTT i. and ( TBI . Note that in the 



Z = 



UL pg 



on 

1 

2A2 



n=o,(j=o 



(16) 



We note that the q 2 dispersion in t pg (Q) means that for a range 
of low T, A pg will vary as T 3/2 . We note that, below T c 
the imaginary contribution in Eq. dT5l > Tq — * faster than 
q 2 as q — > 0. It should be stressed that this approach yields 
the ground state equations and that it represents a physically 
meaningful extension of this ground state to finite T. We em- 
phasize that the approximation in Eq. (0 is not central to the 
physics, but it does greatly simplify the numerical analysis. 



B. Detailed Behavior of the Self Energy 

We have seen that, after analytical continuation, the self 
energy is given by S(k, u>) = S sc (k, ui) + S pg (k, oj), where 



E(k,w) 



w + Ck 
A k , sc 



S pg (k,w) 
A? 



W + ^k w + Ck 



(17) 
(18) 



The BCS-Leggett ground state equations 26 follow, provided 
one makes the approximation contained in Eq. (0. In invok- 
ing this approximation we are in effect ignoring the differ- 
ence between condensed and non-condensed pairs which can- 
not be strictly correct. The simplest correction to E pg (which 
should apply above and below T c ) is to write an improved 
form which most importantly accommodates the fact that the 
coherent Cooper pairs of the condensate are infinitely long 
lived, whereas the incoherent or non-condensed pairs have a 
finite inverse lifetime 7 



E pg (k,o;) 



Al 



+ £k + ij 



E(k,w). 



(19) 
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Here E(k, u>) represents the lifetime associated with channels 
other than the pairing channel and, as is conventional, we pa- 
rameterize E(k, 10) = — iEq. Thus we have 



E(k,w) = 



A? 

k,PS 



6c + il 



- «E 



(20) 



The above equation contains a well known form for E pg . 
It also contains the important addition of E sc . The model 
for E pg was determined in the present context on the ba- 
sis of detailed numerical studie a 30 ' 39 and has been deduced 
independently 40 and widely applied. 31 in the cuprate litera- 
ture. Here the broadening 7^0 and "incoherent" background 
contribution Eo reflect the fact that noncondensed pairs do 
not lead to true off-diagonal long-range order. While we can 
think of 7 as a phenomenological parameter in the spirit of 
the literatur e) 31 ' 41 we stress that there is a microscopic basis for 
considering this broadened BCS for m 37 ' 42 . The precise value 
of 7, and its T-dependence are not particularly important for 
the present purposes, as long as it is non-zero at finite T. By 
contrast E sc is associated with long-lived condensed Cooper 
pairs, and is similar to E pg but without the broadening. It 
is, moreover, often assumed that — iEo ~ — ij, although this 
assumption is not necessary. 



C. Spectral function and Superfluid Density 

The resulting spectral function, based on Eq. ([T9T l and 
Eq. (TTTb is given by 



A(k, 



2A^ k7 (c + e k ) 2 



(^ + £k) 2 (^ 2 -££) 2 +7 2 (^ 2 -£: 



k A sc,k) 2 

(21) 

For convenience, here we do not show the effects of the E 
term. Above T c , Eq. (f2Tb is used with A sc = 0. It can be seen 
that at all k and below T c , this spectral function contains a 
zero at uj = — whereas it has no zero above T c . This means 
that a clear signature of phase coherence is present when one 
passes from above to below T c , as long as 7 ^ distinguishes 
the non-condensed from the condensed pairs. 

These dramatic effects of the condensate in the spectral 
function are also important for addressing the specific heat 
jump at T c which must be present as a thermodynamic indica- 
tion of the phase transition. The onset of a condensate below 
T c (with no lifetime broadening, 7 = ) in contrast to the 
lifetime broadened contribution from the pseudogap is associ- 
ated with clear signatures in the specific heat 30 as the system 
develops superconducting coherence. 

Physically, one can anticipate that the non-condensed pairs 
represent an additional mechanism for destroying the conden- 
sate. It is important to stress that as a consequence this ap- 
proach is different from a Fermi liquid based superconductor 
which has often been presumed in the theoretical literature 9 . 
Because the normal state is, by consensus, a non-Fermi liquid, 
and because there is a smooth evolution from above to below 
T c , it should not appear surprising that the superconducting 



phase is also non-Fermi liquid-based. Important to the analy- 
sis of the superfluid density is the imposition of gauge invari- 
ance through a Ward identity. In this way one find a 36,43 that 
the pseudogap contributions via E pg to the superfluid density 
precisely cancel, in contrast to those from A sc . 

After this cancellation, the superfluid density is found to be 
of the simple form 3 ^ 



n s (T) 



A 2 g (T) 
A 2 (T) 



' n s (T,A(T)) \ 
m J 



BCS 



(22) 



Here, importantly, the quantity (n s (T, A{T))/m) cs cor- 
responds to the conventional BCS form for the <i-wave 
superfluid density, albeit with an unusual, essentially T- 
independent gap A(T) in the underdoped regime. In sum- 
mary, one sees that n s is additionally depressed by bosonic 
fluctuations which insure that n s vanishes at T c , not T*. 



D. Abbreviated Model 

To make the present formalism more widely accessible we 
construct a simplified or abbreviated model in which T* and 
T c are effectively fit to the cuprate phase diagram and the vari- 
ous gap parameters A pg and A sc which appear in the spectral 
function are then readily deduced. For the purposes of the 
present paper we do not focus on this short cut scheme, but it 
serves to make the results here easily reproducible by others. 

We have seen that in the temperature regime below or only 
slightly above T c , the thermodynamical energy gap A(T) and 
its component A pg (T) satisfies A 2 (T) = A 2 g (T) + A 2 C (T) 

where we define E^ d = \J (££ d ) 2 + A 2 , and presume that 

A(T) = A m f(T) satisfies the (two dimensional, mean field) 
BCS gap equation 







k 

\3/2 A 2 



2El d 



(T/T C )^A 2 (T C ), 
A 2 (T), 



T<T C , 
T>T C . (24) 



Here the superscript 2d refers to the fact that we drop the 
third dimension in the energy dispersion so that t z — > 0. At 
each x, the parameter U is chosen to yield the measured T* 
and, knowing T c , A(T C ) can be determined. These equations 
must be solved in conjunction with a self consistent particle 
number equation for [i. Lying behind this phenomenological 
approach is the fact that in a fully consistent theory, 34 . T c is 
(logarithmically) dependent on the inter-layer hopping t z , and 
it vanishes when this parameter is absent where the system is 
strictly two dimensional. Thus we can view t z as a fitting pa- 
rameter which depends on hole concentration x. In the fully 
self consistent scheme one recovers the entire cuprate phase 
diagram for T*(x) and T c {x) by a proper choice of U(x) and 
t z (x). The short cut scheme then allows one to calculate with- 
out too much effort, the various gap parameters as a function 
of temperature and x which appear in the spectral function. 
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Figure 3: (Color online) Spectral function A(<j>, u>) at T/T c — 
1.1, 0.9, 0.1 (from top to bottom) for = 9° (black) and <j> = 36° 
(red). Black and red arrows indicate size of the spectral gap, which 
is measured in ARPES. 
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Figure 5: (Color online) Parameter Insensitivity. This is illustrated 
for T/T c — 0.9. Here we restricted 7 to produce appropriately 
large arcs in the normal phase. Within this range there is virtually 
no change in the size of the deduced spectral gap. We explore two 
orders of magnitude variation in So and again find no change in the 
spectral gap size. 



III. NUMERICAL RESULTS 



A. General Properties of the Spectral Functions 



Figure 4: (Color online) ARPES gap (blue thick line), Ak (red thin 
line), and Ak. sc (green thick line) as a function of <j) at T/T c — 0.9. 



We see that because the total gap A(T, x) satisfies the BCS 
equation there is a BCS-relation between T* and A(T = 0). 
In this way Eq. d23l implies that the excitation gap A contains 
the energy scale T*, not T c . Indeed, at intermediate values 
of the attractive interaction \ U\, A(T) is essentially indepen- 
dent of temperature from the ground state (where A pg = 0) 
to well above T c . We will not discuss the hole concentra- 
tion dependence x in detail in this paper, because it has been 
treated elsewher e) 32 ' 36 . Finally, we note that within this BCS- 
BEC scenario, the mechanism for pairing need not be speci- 
fied. Nevertheless, it is clear that the increase of T* with de- 
creasing x requires that the attractive pairing interaction must 
become stronger as the Mott insulator phase is approached. 



We turn now to detailed numerical calculations of the be- 
havior of the spectral function, A(<p, u>) on the Fermi surface 
(where e-^ — n — 0). Throughout we will define the spec- 
tral (or ARPES) gap as one-half the peak to peak separation 
in the spectral function (when it exists). The dispersion ek 
is obtained using our two dimensional tight binding model. 
For the most part we will consider a prototypical hole con- 
centration x = 0.125, which is associated with a particular 
value of U in Eq. CO} leading to T c /T* « 0.5. We choose 
a bandwidth of At = 250 meV and this results in a T = 
gap about 34 meV. Our results are insensitive to the specific 
parameter set as we will demonstrate below. The only con- 
straint to be imposed from experiment is that there must be 
sizeable Fermi arcs (of order, say, 10° out of 45°) in the nor- 
mal phase, for a moderately underdoped sample. This means 
that the parameter 7 at T c is not much less than about one 
half A at the same temperature. The parameter £0 is found 
to be relatively unimportant for the purposes of the plots we 
present here. It is reasonable to presume that the lifetime of 



7 




Figure 6: Contrasting nodal and anti-nodal temperature dependences in the d-wave case. Figure on the left is the ARPES gap as a function of 
angle </> at T/T c = 1.1, 0.99, 0.1 (labeled on the figure). This figure should be compared with the experimental plots on the right taken from 
Figure 4b in Ref.0 



the non-condensed pairs increases as temperature is lowered, 
since their number becomes fewer. For definiteness, following 
Ref.HH we take Eo = 26 meV independent of T and 7 = 26 
meV at 95 K with 7 (T) = 7(95 K)(T/95 K) above T c and 
7 = 7(T C )(T/T C ) 3 below T c . To be more consistent with 
experimental data, when spectral functions are presented we 
convolve the spectral function with a Gaussian instrumental 
broadening curve with a standard deviation a = 3 meV. 

Figure [3] illustrates the temperature evolution of the spec- 
tral function for <fi — 9° (close to the antinodes) and = 36° 
(close to the nodes) at T/T c = 1.1, 0.9, 0.1 from top to bot- 
tom. Above T c (top panel) the well understood behavior^i^ 
sets the stage for the normal phase which underlies the su- 
perconducting state in the next two panels. In this top panel, 
one sees Fermi arcs, which derive from the broadening term 7 
in E pg , in the near-nodal direction, and a pseudogap in the 
spectral function, associated with A pg near the anti-nodes. 
These arcs appear over that range of k values for which 7 
is larger than the momentum dependent pseudogap. When T 
is slightly below T c (middle panel), a dip in the spectral func- 
tion at = 36° suddenly appears at to = 0. At this <fi the 
underlying normal state is gapless so that the onset of the ad- 
ditional component of the self energy via E sc with long-lived 
pairs (7 = 0) leads to the opening of a spectral gap. 

By contrast, the presence of this order parameter is not re- 
sponsible for the gap near the anti-nodes ((/> = 9°), which, 
instead, mostly derives from A pg . Here the positions of the 
two maxima are relatively unchanged from their counterparts 
in the normal phase. However, A sc does introduce a sharp- 
ening of the spectral function, associated with the deepening 
of the dip at u> = 0. This can be seen analytically from 
Eq. (|2TT > by noting that E sc suppresses A(uj) near u> = 0. 
When T <C T c (lower panel), pairing fluctuations are small so 
that A(T) « A SC (T) and one returns to a conventional BCS- 
like spectral function with well established gaps at all angles 
except at the precise nodes. 

It is useful to look at the behavior of the ARPES gap over 
the entire range of <fi, as studied experimentally- 7 .. To empha- 
size that the spectral gap does not precisely correspond to the 
self energy gap components, in Fig. [4] we plot the spectral 



function gap along with A, and A sc as a function of angle 
at T/T c = 0.9. The figure illustrates that, near the anti- 
nodes, the spectral gap reflects the magnitude of A. Near the 
nodes, however, the spectral gap is more directly associated 
with A sc , in the sense that this gap appears only in the ordered 
phase. The second of these observations is in line with previ- 
ous experimental findings 6,8 . However, it has generally been 
assumed that at the anti-nodes the behavior is governed by the 
so-called "pseudogap". We stress that our interpretation is not 
at odds with this literature. Rather we refer to the full gap at 
the anti-nodes as A(T) which is roughly a constant in temper- 
ature. This contains two contributions, one from A pg (T) and 
one from the order parameter A SC (T). While near T c the for- 
mer dominates, near T w 0, the latter is the more important. 
Thus the gap at the antinodes reflects superconducting order 
as well, at least in these moderately underdoped cuprates. 

Figure[5]shows that the spectral gap shown by the blue lines 
in the previous figure for T/T c = 0.9 is only very slightly 
modified when the parameters Eo ( m the top panel) and 7 (in 
the bottom panel) are altered. While the height of the peaks 
in the spectral function plots will be affected, the important 
derived quantities such as the spectral gap plotted in the figure 
are not changed when Eo is varied by two orders of magni- 
tude. Moreover, if 7 is reasonably constrained to yield a size- 
able Fermi arc above T c , then the behavior of the spectral gap 
below T c does not depend on the detailed values for 7. 



B. Comparison Between Theory and Experiment 

Recently there has been an emphasis on experiments which 
contrast the behavior around the gap nodes with that around 
the gap maxima (or anti-nodes). The right panel of Figure [6] 
indicates the size of the ARPES or spectral gap as deduced 
from one-half of the peak to peak separation in the spectral 
function. These data 7 address a moderately underdoped sam- 
ple. The three different curves correspond to three differ- 
ent temperatures with the legend the same as that in the left 
panel (representing the results of the present theory.) Impor- 
tantly, one sees a pronounced temperature dependence in the 
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Figure 7: (Color online) (a) The ARPES gap (red squares), Ak (thick 
blue dashed line), Ak, sc (black dashed line), and Ak, P9 (orange dot- 
dash line) as a function of T/T c for <j> = 36°. (b) The ARPES gap 
as a function of T/T c for </> = 36° (red squares) and <f> = 30° (green 
circles). This panel should be compared with Figure 2d of Ref.tTJ. 



behavior of the ARPES spectral gap for the nodal region (near 
45°), as compared with the anti-nodal region (near and 90°), 
where there is virtually no T dependence. 

Theory (on the left) and experiment (on the right) are in 
reasonable agreement and one can readily understand the con- 
trasting temperature response associated with the different k 
points on the Fermi surface. To see this, note that the nodal 
regions reflect extended gapless states or Fermi arcs 24 above 
T c . It is natural to expect that they are sensitive to the onset of 
A sc , in the same way that a strict BCS superconductor, (which 
necessarily has a gapless normal state), is acutely sensitive to 
the presence of order. By contrast, the anti-nodal points are 
not as affected by passing through T c because they already 
possess a substantial pairing gap in the normal phase. 

The dramatic variation in the temperature dependence of 
the spectral gap as one moves along the Fermi surface has 
given rise to the so-called "two gap scenario" 1 . In (perhaps) 
overly simplistic terms the one gap and two gap scenarios 
are differentiated by the presumption that in the former the 
pseudogap represents a precursor to superconductivity, while 
in the latter the mysterious cuprate pseudogap is viewed as 
arising from a competing order parameter. The two gap sce- 
nario is viewed as a consequence of a number of different 
experiments 113 all of which have been interpreted to suggest 
that the antinodal region is associated with this alternative 
(hidden) order parameter pseudogap and the nodal region is 
dominated by superconductivity. By contrast the viewpoint 
expressed here (based on BCS-BEC crossover theory) leads 



naturally to a different T dependence for the nodal and anti- 
nodal region, but at the same time it belongs to the class of the- 
ories which argue that the pseudogap is intimately connected 
with the superconductivity. 

We turn in Fig. [7] to very important temperature dependent 
studies 7 which suggest that the nodal gap may directly reflect 
the order parameter. Figure [3 a) plots the various gap param- 
eters, A(T), A P9 (T) and A SC (T) in the self energy as com- 
pared with the spectral gap measured near the node at <fr = 36° 
(indicated by squares) as a function of temperature. It can be 
seen that this spectral gap, while it is distinct from the order 
parameter A SC (T) (except at the lowest temperatures), van- 
ishes rather close to T c . The figure shows that the gap parame- 
ter A(T) is relatively constant through T c , so that the decrease 
in A pg (T) with decreasing T is compensated by the increase 
in A SC (T) through the inter-conversion of non-condensed and 
condensed pairs. To compare directly with experiment, in 
Fig- 13b) we plot the spectral gap for two different angles, cj>, 
as a function of T, in a fashion which looks rather similar to 
Fig. 2(d) of Ref. 1, For tj> = 30°, which is somewhat further 
from the nodes there is a small spectral gap (pseudogap) above 
T c . Because of the tpk factor, closer to the anti-nodes the over- 
all magnitude of the ARPES gap is larger than at cj> = 36°. 

In Fig. [8] we address the important issues which have been 
raised in Refs. 0, 0. andl24l These papers make the case that 
the pseudogap is a consequence of the superconductivity. The 
figure in the main body is a plot of the spectral gap for a few 
different temperatures from above to below T c as a function 
of the simplest <i-wave form for tp^. This figure compares 
favorably with Fig. 3(b) in Ref. 7. The central point illustrated 
here is that at the lowest temperatures one reverts, in effect, to 
a simple one-gap scenario. That is, the BCS-like ground state 
wavefunction obtains with A = A sc . 

In the inset of Fig. [8]we present a contour plot of the occu- 
pied spectral weight corresponding to the product of the spec- 
tral function and Fermi function. In this way one can infer 
the dispersion relationship associated with the normal phase 
and see to what extent it is related to that below T c . The left 
panel is below T c and the right panel above T c . This contour 
plot, albeit represented differently, compares rather favorably 
with Fig. 4 in Ref. 2. The similarity of the two panels would 
not be expected if the pseudogap were related to another order 
parameter. 

Together Figure [8] and related experiment s 2 ' 7 ' 24 provide ev- 
idence that the pseudogap has to be viewed as ultimately as- 
sociated with the superconductivity. The normal state excita- 
tions appear to have a (broadened) BCS-like dispersion. The 
nodal and anti-nodal behavior appear to be intimately con- 
nected in the ground state. 



IV. PHENOMENOLOGICAL MODEL FOR HEAVILY 
UNDERDOPED SYSTEM 

There is a growing body of work on more heavily under- 
doped cuprates 6,8,13 from which one can infer that the simple 
d-wave, BCS-like ground state may not be appropriate nearer 
to the insulating phase. Here, if one looks at the experimen- 
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Figure 8 : The ARPES gap as a function of | cos ( k x ) — cos ( k y ) | /2 for 
T/T c = 0.1 (blue solid line), 0.99 (green dashed line), and 1.1 (red 
dot-dash line). This should be compared with Figure 3b of Ref. 7. 
Inset is a contour plot of the occupied spectral weight at <f> = 22.5°, 
showing peak sharpening below T c . We follow a similar sweep as 
that in Ref. |2J and the white line indicates the intersection with the 
Fermi surface. Here the intensity corresponding to below (left panel) 
and above (right) T c is largest( smallest) in the red(blue) and we have 
taken smaller 7 for illustrative purposes. 

tal analogue of Figure[8] the lowest temperature behavior still 
exhibits a deviation from the simple cos k x — cos k y form. In- 
deed kinks are often seen 8 somewhat like that shown in Figure 
[SJ but for the case of very low temperatures. The kinks are as- 
sociated with the fact that the ARPES gap curves in the nodal 
region seem to reflect the superconducting order while, as be- 
fore the antinodal behavior reflects what is referred to as the 
pseudogap. As a result it has been argued that 8 "the very dif- 
ferent properties of these two gaps lead us to conclude that 
there is no direct relationship between the pseudogap and the 
superconducting gap". 

Because there appears to be a rather continuous 7 evolution 
from moderate to heavy underdoping, we, instead speculate 
that the physics of the pseudogap in the two regimes must 
be rather similar and that the non-simple d-wave ARPES gap 
behavior at the lowest temperatures in heavily underdoped 
cuprates is a natural extension of the higher T <T C behavior 
seen at moderate underdoping. At these higher T < T c there 
are two gap components A pg and A sc 7^ 0. Thus, a rea- 
sonable precursor- superconductivity- based phenomenologi- 
cal model for this extreme underdoped regime is to presume 
that A pg persists into the ground state, perhaps because of a 
contamination from the near-by insulating phase. We view 
this insulating state as introducing a finite value for the zero 
temperature pseudogap. This is consistent with the way^ 4 - the 
insulating phase appears in our calculations where the pg gap 
component persists to the lowest temperatures while A sc is 
strictly zero beyond a critical value for the attractive interac- 
tion, or equivalently a critical value for T*. 

We emphasize that all previous discussions and figures have 
been microscopically based and derived, but in this section we 



proceed purely phenomenologically. The goal of this discus- 
sion is to arrive at a model for the extreme underdoped case 
which is smoothly connected to the physical picture we have 
thus far exploited for more moderately doped cuprates. We 
need to incorporate a clear deviation from the d-wave ground 
state, kinks or other breaks in the ARPES gap function which 
distinguish different gap shapes around the nodal and anti- 
nodal regimes, and clear evidence for incoherence even below 
T c , but only near the antinodes. The model we present grew 
out of a discussion with A. Yazdani and his collaborators 44 
who have observed a similar gap shape in their STM experi- 
ments. 

To describe this class of models we assume that all gap 
functions (but not the spectral gaps themselves) have the form 
\ g ,k = A pg ip^ and A SCjk = A sc ^ k . At a given tempera- 
ture, the pseudogap now has two contributions: one from the 
usual preformed pairs, which will ultimately go into the con- 
densate at sufficiently low T and another from the admixture 
of insulating state which we view as a "zero temperature pseu- 
dogap". In this way there is a weak temperature dependence 
in A pg associated with the pair conversion process and con- 
comitantly A sc is also T dependent. A typical parameter set 
is shown in the inset of Figure [9] This plot is to be contrasted 
with the behavior shown in Figure [7^. 

For definiteness we presume that the total excitation gap 
is given by the mean field gap A m f (T) defined in Eq. d23l ), 
so that the superconducting order parameter contribution is 

A SC (T) = ^jAl lf {T)-A 2 pg (T). The pseudogap contribu- 
tion is written as A pg {T) = yj A^ g0 (T) + A^ gl (T) with the 
zero temperature pseudogap given by A pg o (T) = aA m f (T) 
and A Ml (T) = (T/Tc) 3 / 2 Jtf~(T) - A 2 g0 (T) for T < 

T c and A pg (T) = A mf (T) for T > T c . Here A m/ (T) 
is the gap obtained from a mean-field calculation of d-wave 
BCS theory, as derived from Eq. (f23t . In our microscopic 
calculations one would have a = which appears consistent 
with moderately underdoped systems. However, for heavily 
underdoped cuprates We choose a such that the sc and pg 
contributions at T=0 are in the ratio of 1 :2, as a typical exam- 
ple. We take 7 (T) = E (T) = 0.5A p9 (T) for T < T c and 
7 (T) - E (T) - (T/T C ) 7 (T C ) forT > T c . 

Figure [9] shows the behavior of the spectral gap for the 
heavily underdoped model and for a range of temperatures 
T/Tc = 0.1, 0.8, 0.9, 0.99, 1.1. The dashed line is an extrap- 
olation of the simple ci-wave fitted form found near the gap 
nodes and associated with the order parameter A sc at the low- 
est temperature. While there is a simple d-wave fitted form 
also at the anti-nodes the effective gap here is the much larger 
parameter A, which consists mostly of a pseudogap contribu- 
tion, for this heavily underdoped system. 

Figure [10] shows a plot of the actual spectral functions at 
two angles <\> — 36° in red and tf> = 9° in black at three dif- 
ferent temperatures from above to just below T c to finally at 
T/T c = 0.1. The behavior in this heavily underdoped system 
can be contrasted with that shown in Figure [3] for moderate 
doping. The nodal curves show the Fermi arc behavior above 
T c , followed by the opening of a gap (which reflects super- 
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Figure 9: (Color online) Behavior of the spectral gap as a func- 
tion of angle <j> for a phenomenological model representing a heav- 
ily underdoped system. The inset plots the gap functions which 
should be contrasted with that shown in Figure 7a. The dashed 
line is the extrapolation of the simple cos(2</>) behavior found near 
the gap nodes. Solid curves from top to bottom correspond to 
T/Tc = 0.1, 0.8, 0.9, 0.99, 1.1. 
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conductivity) below T c and the ultimate establishment of well 
defined coherence with decreased T as evident by the narrow, 
well defined peaks. By contrast the anti-nodal regime (unlike 
its counterpart in Figure [3) does not indicate the presence of 
coherent quasi-particles. Rather, even at the lowest tempera- 
tures the peaks are broad, and very little changed from those 
above T c . 

There are features of this model which do not capture all 
the phenomena observed experimentally. The "kink" effects 
seem to be strictly associated with the Fermi arcs of the nor- 
mal state and not particularly close to the magnetic zone 
boundary 6 , since the arc size is rather small in this underdoped 
regime. Moreover we have presumed a strictly e?-wave gap 
shape which constrains the behavior of the spectral gap near 
the antinodes. Nevertheless, this is a reasonable model for fur- 
ther study, since it does preserve some of the key physics of 
the experiments. 



V. CONCLUSIONS AND COMPARISONS WITH THE 
LITERATURE 

This paper addresses issues which are at the center of 
major debates in high temperature superconductivity. Do 
the recent (so-called "two gap" ) experiments which report 
a difference associated with the nodal and anti-nodal re- 
sponse in ARPES&iS or in Raman^S or scanning tunneling 
microscopy 1314 rule out the possibility that the pseudogap de- 
rives from the superconductivity itself? We argue that, despite 
strong claims in the literature, pseudogap formation owing to 
preformed pairs is, in fact, consistent with these experiments. 
We stress that our approach for the moderately underdoped 
cuprates is not phenomenological. It was in place well 30 be- 
fore these experiments were undertaken. 



Figure 10: (Color online) Spectral functions with convolution for 
phenomenological model of a heavily underdoped system. This 
model should show that at the lowest T the behavior around the 
antinodes is not much more coherent than that in the normal state. 
This figure should be contrasted with Figure[5] 



We have emphasized that our explanation for the physics 
is relatively simple and is based on a stronger-than-BCS at- 
tractive interaction, associated with short coherence length 
Cooper pairs. The formation of isolated pairs (in contrast 
to extended regions of fixed pairing amplitude) takes place at 
T*, while condensation appears at T c . What is crucial is that 
pseudogap effects which are associated with these pre-formed 
pairs do not disappear immediately below T c . Rather they 
persist as non-condensed pair excitations of the condensate. 
This is not a Fermi liquid based form of superconductivity, 
because there are bosonic degrees of freedom associated with 
the fermion pairs. Nor should this be thought of as a "one 
gap" picture. There are two components to the pairing gap, 
one from the non-condensed pairs and another from the con- 
densate. 

A central equation is Eq. ( |20b which shows that both com- 
ponents are important in the self energy and therefore in the 
spectral function. The contribution from the preformed pairs 
T, P g is crucial for forming the Fermi arcs above T c . These ap- 
pear in the nodal regions where 7 is relatively larger than the 
momentum dependent gap. The contribution from the con- 
densate E sc is crucial just below T c because it opens up a true 
gap in the Fermi arc region. This is reminiscent of a con- 
ventional BCS superconductor which necessarily has a gap- 
less normal state and is, thus, extremely sensitive to the pres- 
ence of coherent order. This is, in contrast to the anti-nodal 
regimes where the large pseudogap above T c is very little af- 
fected by the addition of the superconducting order, except 
through peak sharpening or coherence effects. 
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In the context of Eq. ( f20b it is generally believed 31 that there 
is only one component to the self energy (£ M ) and that the 
onset of coherence coincides with a dramatic decrease in 7 
below T c . We strongly disagree with this assumption. Rather 
there are two contributions to the self energy below T c and 
only one above. Thus, one should not argue that 7 precisely 
vanishes at T c but rather there is a continuous conversion from 
non-condensed to condensed pairs as T is lowered within the 
superfluid phase. The non-condensed pairs below T c have fi- 
nite lifetime while the condensed pairs do not. 

In this paper we also discussed the fact that there is support 
from another class of experiments that the pseudogap and the 
superconducting gap are intimately connecte d 2 ' 7 ' 24 . The low- 
est temperature spectral properties 7 — of, at least, moderately 
underdoped samples seem to fit a simple d-wave angular de- 
pendence and recent normal state data 2 provide evidence for a 
dispersion deduced from the spectral function which is similar 
to that in the superfluid phase. 

Finally, we addressed heavily underdoped cuprates in a 



phenomenological fashion. Here the simple d-wave gap shape 
may not be appropriate 8 . We argued that what is crucial is 
that there is a continuous evolution from moderate to extreme 
underdoping 7 so that it is unlikely that the pseudogap has a 
different origin in the two regimes. Rather some of the same 
physics must be at play. We postulated that there may be 
a zero temperature pseudogap present in highly underdoped 
systems which may derive from some admixture of the insu- 
lating phase. 

In summary, this paper has shown how to reconcile a wide 
class of experiments in the moderately underdoped cuprates 
within a pre-formed pair framework where there are, neverthe- 
less, two components to the energy gap. This framewor k 33,36 
predates the class of experiments we address here. 
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Davis, A. Yazdani, Colin Parker and Aakash Pushp, as well 
as Wei-Sheng Lee and D. Morr for helpful discussions. 
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